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Cylindrical Antenna Theory l 
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A partial survey of cylindrical antenna theory pertaining to a tubular model with a 
narrow gap is presented. The survey includes discussion of the theories of Hallen, King and 
Middleton, Storm, and Zuhrt. A conceptual relation between theory and experiment is 
described. The latter part of the article is concerned with a new Fourier series solution 
of the Hallen equation. This solution is developed in such a way that the expansion coeffi- 
cients are the unknowns of a system of linear equations. The elements of the coefficient 
matrix are given by a highly convergent series. Numerical results are given for half and 
full wavelength antennas with half length to radius ratios of 60 and 5007T. These results 
compare quite closely with those obtained from King-Middleton theory. 



1. Introduction 

Existing solutions to Hallen's integral equation 
for the current distribution on cylindrical antennas 
fall into two main categories [1, 2]: 2 

1. Iterative solutions which use an approxima- 
tion to cylindrical antenna current as a starting 
point. Successive iterations generate improvements 
in the original assumption. Approximations are re- 
quired at some stage in the process if tractable inte- 
grals are to be obtained. The approximations are 
not severe if h/a, the ratio of antenna half length to 
radius, is large. Impedances obtained from itera- 
tive solutions are in good agreement with experi- 
ments performed on thin antennas. Although suc- 
cessful, iterative solutions become laborious beyond 
second or third order, and the approximations 
become suspect in the case of thick structures. 

2. Solutions in which the integral equation is con- 
verted into a set of linear simultaneous equations 
with Fourier coefficients of the current distribution 
as unknowns. Typical of these are the theories of 
Storm and Zuhrt. (Strictly speaking, Zuhrt did 
not solve Hallen's integral equation, but one deriv- 
able from a somewhat different point of view.) 
Storm approximated matrix elements in his set of 
equations so that his theory is limited to thin struc- 
tures. In addition, Storm's theory contains two 
rather fundamental errors which, acting in concert, 
produce fortuitous results. Zuhrt obtained matrix 
elements by graphical integration, a sufficiently tedi- 
ous process to limit his calculations to low order. 
Neither of these solutions fully exploits the Fourier 
series technique for obtaining the current distribut ion. 

We have obtained a solution similar to Storm's in 
which computation of matrix elements can be easily 
done with high accuracy even for h/a ratios as low 
as 8 or 10. Matrix inversion is easily accomplished 
with modern digital computers so that solutions of 
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high order are feasible. We have carried out cal- 
culations to 25th order for half and full wavelength 
antennas with h/a ratios of 60 and 500Y. Results 
from the King-Middleton iterative solution com- 
pare favorably to ours so that their work for h/a > 
60 has been checked by comparison with an exact 
theory. 

Although this paper is principally concerned with 
developing a Fourier series solution, a reasonably 
complete treatment of iterative solutions is included 
in an attempt to provide a self-contained account 
of cylindrical antenna theory. Even with this aim 
in mind, the treatment of iterative solutions, as 
well as the theories of Storm and Zuhrt, is sufficiently 
involved that the reader is referred to the original 
work for many' of the details. 

2. Statement of the Problem 

Vector potential as a function of position, A(f), 
is given in terms of a current distribution, J(r) } by 
the expression 

The vector integration in (2.1) must be taken over 
all sources of the vector potential field and it is this 
requirement that causes a direct solution of a realistic 
antenna problem to be exceedingly difficult. In the 
case of a cylindrical antenna the integration would 
have to include currents in the antenna, the feeding 
transmission line, and the oscillator which supplies 
power to the antenna-transmission line assembly. 
Every possible combination of antenna, transmission 
line and driving generator would have to be treated 
as a special case and a mathematical solution of any 
given case would by itself be formidable. An ideal- 
ized problem can he extracted from this situation by 
consideration of an extremely thin walled tube of 
infinite conductivity with a narrow circumferential 
gap corresponding to the antenna terminal zone. 
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Choosing the conventional cylindrical coordinates, 
(p, 4>,z), the antenna is defined by p=a and |s| <h. 

A purely hypothetical generator is assumed such 
that the electric field in the gap is azimuthally 
symmetric. The voltage across the gap is defined by 



c/ga 



E s dz, 



(2.2) 



where E s is the ^-component of electric field at p=a. 
E s is zero outside the gap since the simplifying 
assumption of perfectly conducting tube walls has 
been made. If the gap width is decreased as V is 
held constant, we must express E s by 



E 8 =-V8(z) for \z\<h, 



(2.3) 



where 5 (z) is the usual delta function. Center-fed 
models will be considered here although this restric- 
tion can be removed. 

The value of electric field anywhere in space is 
given by 

E=— grad 3>— juA, (2.4) 

where $ is the scalar potential. (Time dependence 
proportional to e jut is assumed in all of the relation- 
ships used here.) Scalar potential can be found 
from the Lorentz condition 



<£>=— div A. 

cofxe 



(2.5) 



The combination of a tubular model without end 
caps and the assumption of a symmetric field leads 
to the conclusion that only the 2-component of A 
is different from zero since all of the current sources 
on the tubular surface will be in the ^-direction. 
Equations (2.4) and (2.5) can be combined to give 



E 2 = — J — (c> 2 A/c>z 2 +k 2 A), 
cojue 



(2.6) 



where the symbol A without the vector bar simply 
stands for A z , and P = w 2 pe. In general A is a 
function of both z and p. If the operator (5 2 /ds 2 +P) 
is applied to A (z, p) and then p is taken at the 
antenna surface, we obtain 



Cx)(JL€ 



b 2 A s /dz 2 +k 2 A s =^V8(z), 



(2.7) 



an equation for the surface value of vector potential 
valid for \z\<h. Solutions of the homogeneous 
equation, d 2 A s /dz 2 + k 2 ^4 S =0, are simply cos kz } sin 
kz, e jicz or e~ jkz . Linear combinations of these solu- 
tions may be used to build a solution to (2.7). Thus, 



or 



A s =d cos kz+Di sin k\z\, (2.8) 

A S =C 2 cos kz+Dsfi-M'K (2.9) 

These solutions are completely equivalent. How- 



ever, we prefer (2.9) as a basis for studying certain 
properties of the infinite cylinder. D x and D 2 are 
evaluated by substituting these solutions into (2.7) 
and performing the indicated operations with due 
account being taken of the discontinuous derivatives 
of sin k\z\ and e~ jJc]z K The results are D l =(co /j,eV)/2jk 
and D 2 = (aneV)/2lc. Then 

A S =C 1 coskz+^~sm k\z\, (2.10) 



jk 2 
A 8 =C 2 cos kz+^f- ^e-M*K 



(2.11) 



Since e~ jkl ^ can be written as cos kz— j sin k\z\, it is 
easily shown that Ci—C 2 +o)neV/2k. Either expression 
may be used, the choice between them being only a 
matter of taste. 

Under the assumptions being made, the current 
distribution of the general formula (2.1) degenerates 
to a surface distribution, K(z ). It must be borne 
in mind that the tube has both inner and outer 
surfaces, and K(z ) is the sum of current densities 
on both surfaces. The field point, r, can be taken 
at the surface of the cylinder p=a so that an alternate 
formula for A s is 

A,=£\ K(z ) e adz d<t> , (2.12) 

47rJ_ A J_ 7r \r s —r \ 

where r ranges over the antenna surface during the 
course of the integration. 

Several changes in notation are convenient at this 
point. The field point is (a, z, <f>) and the source 
point has the coordinates (a, z , <t> ). Only the 
difference <f>— </> is of any significance because of the 
azimuthal symmetry. Therefore, we may set <j> to 
zero, change <t> Q to and indicate the angular inte- 
gration as already performed. The quantity z is 
changed to f so that the subscripts may be avoided 
in future formulas. Total current is given by 
/ (f) = 27r aK (f). The quantity \r s — r Q \ is replaced 
simply by R. With these changes 



where 
and 



A * = tS-l m9{z ~ m ' (2 - i3) 



1 C v e~ jkR 
J? = [4a 2 sin 2 0/2+(3-f) 2 ] 1/2 . 



(2.14) 
(2.15) 



We now have two formulas for the surface value 
of vector potential. One has been derived from the 
generic formula for vector potential in terms of cur- 
rent sources by specializing to the geometry of the 
problem under consideration. As such it is a 
general formula for vector potential generated by a 
(^-symmetric current in the 2-direction on a tubular 
conductor no matter what other conditions are to 
be imposed on the problem. On the other hand, 
(2.10) provides a vector potential which leads to the 
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boundary values of electric field desired in the 
present problem. If these expressions are equated, 
an integral equation for the current distribution as- 
sociated with the chosen boundary conditions results: 



4w 



J "l^giz-^di^C, cos kz+^ \ sin k\z\- 

(2.16) 



It is convenient to multiply the previous equation 
through by jk/ujjie and make the following definitions: 



and 



cojue cojue 47r 47T 47r 

§/(r)=/(r). 



No generality is lost by letting V=l. This com- 
pletes the mathematical formulation of the problem. 
We are to consider the solution of 






r)ff(a-rW=C F coBt»+gainifc|»| J |a|<k. 



(2.17) 



The constant C must be determined by the boundary 
condition f (±h) = 0. 

There is no doubt that the integral equation cor- 
responds exactly to the chosen model. It has already 
been pointed out that the model does not correspond 
to any physically realizable antenna. Physical 
antennas may be either solid or tubular conductors, 
and they may or may not be fed in such a way as to 
preserve ^-symmetry. Lack of symmetry (as 
exemplified by a linear antenna fed by a two wire line) 
can be rationalized to some extent if ka is small. 
However, the most serious point is the highly ideal- 
ized nature of the generator region of the mathe- 
matical model. The infinitesimal gap is really a short 
circuit across which a hypothetical but finite voltage 
has been impressed. Thus, the input current and 
admittance of the model are certainly infinite. 

It is not immediately clear that an infinite ad- 
mittance model can be managed mathematically in 
such a way as to yield a physically significant finite 
result. Wu and King have discussed this point in a 
recent paper [3]. They have shown that the singu- 
larity in I(z) near z=0 is logarithmic and of very 
short range. Thus, according to Wu and King, 
"since . . . the singularity actually gives a contri- 
bution to the current distribution only in an exceed- 
ingly small and physically meaningless distance of the 
order of magnitude h exp (—1/ka), it may, in prinei pie 
simply be subtracted out." According to this line of 
reasoning iterative solutions of the integral equation 
are successful because they are started with a con- 
tinuous approximating function and are carried to 
such low order that the singularity does not develop. 

An important aspect of any theory is its relation- 
ship to experiment. It is possible to avoid the in- 
herent singularity in the current distribution which 



is a solution to Hallen's integral equation and obtain 
a finite answer for the theoretical input admittance 
to a cylindrical antenna. The finiteness of the result 
does not, however, guarantee that it is physically 
significant. Experimental antennas must be pro- 
vided with a realistic terminal zone which is con- 
nected to a transmission line of some sort. Measured 
impedances are then complex combinations of 
antenna and terminal zone effects. One way of 
extricating these effects is to make a theoretical 
correction for terminal zone effects on a sequence of 
experimental impedances and extrapolate the cor- 
rected data to the limit of a small terminal zone. 
The residual impedance is then supposed to be char- 
acteristic of the antenna itself and it is this idealized 
inference from experiment rather than raw data 
which is to be compared with theory. 

An outline of a feasible experimental program will 
be helpful at this point. Consider the experimental 
arrangement shown in figure 1. Symmetry allows 
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Figcre 1. Experimental model of a cylindrical antenna. 

us to place a large ground plane at 2=0 and simulate 
half of a center fed antenna. Measuring instru- 
ments are located below the ground plane so that the 
antenna is shielded from extraneous interaction with 
the apparatus. Appropriate measurements are made 
so that the experimenter can determine Z T) the 
equivalent impedance terminating the coaxial line. 
It should be possible to relate Z T to parameters of 
the radial transmission line and to Z a (w), the input 
impedance of a cylindrical antenna with gap width 
2 w. If so, one should be able to calculate Z a (w) 
given measurements of Z T . Z a (w) can be extra- 
polated to obtain Z a (0). Z a (0) obtained in this way 
is independent of the terminal configuration actually 
used. We will not attempt to give a mathematical 
treatment of the experimental situation. It is suffi- 
cient for our purposes to establish a conceptual 
foundation for relating theory to experiment as a 
prelude to subsequent theoretical development. 

Consider now the sense in which a Fourier series 
solution for I(z) "subtracts the singularity. " Since 



571 



the unknown function in (2.17) is even, it is developed 
to finite order as 



/ W (f)=§+Z) F n cos n^/h. 

* 71=1 



(2-18) 



A procedure for determining the F n will be exhibited 
later. At this point it is sufficient to recall that the 
sum of a finite number of terms of a Fourier series 
fits the function being described in a least squares 
sense. Specifically, 

<Z 2 =J J/*(f)-/(f)Nf (2.19) 



is a minimum when the coefficients of /#■(£) are 
Fourier coefficients. Intense, short range variations 
cannot contribute appreciably to f N {$) unless N is 
made so large that even the singular behavior of 
/(f) begins to develop. 

A numerical estimate of the range of the singular- 
ity is worthwhile at this point. Consider an antenna 
for which kh=w/2 and h=60a. Then ka^. 025 and 
l/ka^40. Substitution of these numbers into the 
range estimate of Wu and King leads to the conclu- 
sion that a continuous function can fit the current 
distribution except in a small distance equal to about 
he~*°. Even if the thickness is increased until h= 10a 
the singularity is important only over a range of 
about he-*- 6 ~.0015h. Thus it can be seen that /(0) 
as given by 



2 71=1 



/(0)=^+S F n 



(2.20) 



will be apparently well-behaved in calculations of 
practically feasible order even though theoretical 
considerations indicate that the infinite series must 
diverge. 

An extraneous feature of the theoretical model 
becomes apparent when it is compared to a proposed 
experimental model. Current on the inner surface 
of the theoretical model near the feed point is re- 
placed by current associated with a realizable ter- 
minal zone in the experimental model. If ka is 
less than 2.61a, tube modes are below cutoff and 
are rapidly damped out [4]. Formulas for removing 
the inner current near 2=0 from the total current, 
which is a solution to the integral equation, will be 
presented later. As might be expected, the correc- 
tion makes a small difference in theoretical antenna 
admittance. 

Once the singularity near 2=0 and the tube cur- 
rent are understood and properly removed from the 
theory the remaining possibilities for refinement are 
somewhat limited and consist, for the most part, of 
removing the assumption of an infinitesimally thin, 
perfectly conducting tube. Such considerations 
would be supererogatory in view of the many suc- 
cesses of infinite conductivity models in electromag- 
netic theory. 



3. Solution for the Infinite Cylinder 

In subsequent work the current distribution on a 
finite antenna is expanded in a Fourier series. Such 
an expansion cannot be valid for a function contain- 
ing a singularity unless the singularity is integrable. 
One expects that I(z) is singular in the neighborhood 
of the delta generator and that the nature of the 
singularity is independent of antenna length. It 
can be imagined that the current on a finite antenna 
is composed of waves emanating from the feedpoint 
and waves reflected from the ends of the antenna. 
Only the outgoing waves are expected to be singular 
at 2=0 since they are directly associated with the 
delta generator. It is instructive to consider the 
case of the infinite cylinder before proceeding with 
the solution for a finite antenna. This case can be 
solved exactly in integral form and the Fourier 
transform of the current identified from the solution. 
Asymptotic behavior of the transform gives a clue 
to the nature of I(z) near 2=0 even for the finite 
antenna. 

Consider now the integral equation formed by 
equating (2.11) and (2.13), 



£ J^ /(r)jr(s-fl#=Ci cos kz+^f 



V 

k 2 ' 



-#i*i 



(3.1) 



If h becomes infinite there is no mechanism for 
the formation of standing waves on the cylinder. 
In that event C 2 may be set to zero. The integral 
equation becomes 



Zo 
2tt 



J+co 
- CO 



z-r)df=JV*W;|z|<co. (3.2) 



Solution of this equation is easy if we are armed 
with the identities 

0(s-f) = -| f°° JoWHPiMe'^'-Vda, (3.3) 



and 



-* w =f n5 — r**- ( 3 - 4 ) 



The parameter j} is given by 

/3=(P-a 2 ) 1/2 . 



(3.5) 



Each integral is taken in the complex plane of a 
along the real axis from — oo to + °° with a down- 
ward identation at a=—k and an upward indenta- 
tion at a=-\-k. 

The Fourier transform of I(z) is 

/(«)=i f + "i(r)e-**<*r. (3.6) 

■\ Zt J -co 

This definition and the identities for g(z—$) and 
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e~ iJclzl lead to 



- co -y ZtZjq J — co P 



If (3.7) is to be true for all z, 

AkV 1 



/(«)=- 



Tlien the solution of (3.2) is 



r 2&Ff H 

I(S)= ^J- 



9.7'aZ 



(?J ((3a)m*>((}a) 



da 



(3.7) 



(3.8) 



(3.9) 



Since g ;a2 is the only factor in the integrand contain- 
ing an odd part, (3.9) becomes 



7T A) J 



cos az 



We have defined ft as follows: 
P 2 = k 2 -a 2 , 
P=\P\ for «<& and 
0=-il0| fora>fc. 

Therefore, for large a 

P= —ja. 



(la- (3.10) 



(3.11) 



(3.12) 



Certain identities involving cylindrical functions are 
required: 



J ( —jaa) = I ( — aa) = I Q (ad) , 

J* 



(3.13) 



Io(aa)- 



and 



(27raa) 1/2 
K (aa)->(Tl2aa) 1/2 e-<*. 
For the definition of cosine transforms we take 



/.(■ 



»>=v!l 



7(s) cos a- 26?2. (3.14) 



Equations (3.10) through (3.14) imply that, for 
large a 



Ic((*)->2jl— 

Z a 



(3.15) 



Inspection of tables of cosine transforms [5] reveals 
that when a function behaves as \\\z for small 2, its 
cosine transform behaves as -~yJir/2 a- 1 for large a. 
It follows that for small enough z 



I(z)c 



AkaV, 7 
-j 7 In kz- 



(3.16) 



In the above expression k has been selected as a 
multiplier for \z\ to make the argument of the loga- 
rithmic function a pure number. The asymptotic 
behavior of I c (a) is of no help in deciding whether or 
not k is a proper choice for this parameter. However, 
k seems attractive since it is given by k=2ir/\, and 
the wavelength is a natural unit of length in radia- 
tion problems. 

There is also a somewhat obscure reason for choos- 
ing k as a parameter to convert \z\ to nondimensional 
form. Consider another version of (3.15) in which 
the parameter k is retained as a is allowed to become 
large. As soon as a becomes greater than k, T c (a) 
becomes a pure imaginary, J {) ((3a) and H& 2) (pa) go 
over to I 9 (\P\a) and ( — 2/jw) K (\P\a). Then using 
asymptotic forms, 



I(a>k)- 



.2^2ivkaV 



1 



(c?-k 2 ) l/ 



(3.17) 



I c (a) can now be regarded as having been separated 
into two parts: t c (a<^k) which is zero if a^-k and 
I c (a^>k) which is zero if a<^k. Only I c (cT>k) is of 
interest at present since it is responsible for the singu- 
larity in l(z). Now I c (a^>k) is asymptotically pro- 
portional to the transform of the Neumann function, 
Y (kz). The latter contains a singular part which 
is proportional to In k z. 

If the above arguments in favor of choosing (3.16) 
for the form of the singularity are acceptable, we 
can proceed with an estimate of the range over 
which this expression is a good approximation to 
antenna current. For this purpose we require an 
estimate of l(z) outside the range in which the 
logarithmic function is dominant. It is difficult to 
obtain such an estimate from (3.10). However, an 
alternative argument can be constructed. It is 
known that an expression of the form I m sin k(h—\z\) 
is a fair approximation to current on a finite antenna 
over most of its length. It is convenient, rather 
than necessary, to use knowledge of the finite length 
antenna in estimating the range of the singularity. 
We now arbitrarily establish the criterion that the 
logarithmic function is to be used in the range 
H<C*7 where rj is given by 



4kaV 



In kr]^\I m 



(3.18) 



\I m \ is of the order of 0.01 when V=l. Z is of 
the order of 400. An order of magnitude estimate 
for 7] is 



„~-e 



-Ilka 



(3,19) 



This is a small fraction of a wavelength even for 
quite thick antennas. Moderately large changes in 
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\I m \ do not affect this conclusion appreciably. 

Further examination of this question would re- 
quire analytic or numerical inversion of (3.10). 
Either approach is apparently formidable. 

4. Iterative Solutions for Finite Antennas 

A brief discussion of iterative solutions is provided 
here as background for the general reader. We are 
primarily interested in presenting a critique of the 
approach rather than compiling a report on the 
voluminous literature of the cylindrical antenna 
problem. Consequently, the presentation omits 
many points which are essential to a detailed under- 
standing of the iterative method. It does present a 
few features which require comment by way of 
justifying additional consideration of finite antenna 
theory. 

The first task is to cast eq (2.17) into a form 
suitable for iteration. To that end the quantity 

/ 00*00= f* j(z)w{z } S)dt (4.1) 

J -h 

is added and subtracted to the left hand side of 
(2.17). It is convenient to abbreviate by letting 

O cos kz+^smk\z\=P(z). (4.2) 

With these changes and some elementary transpo- 
sitions, (2.17) can be written 

/(2)= ^{ P(2) ~Jj /f) ^ (s_f)_/(2)w(3 ' f)Mf }- 

(4.3) 

Assuming that w(z, f ) has been selected, a program 
for obtaining a solution is: 

1 . Substitute an approximation to j(z) under the 
integral on the right hand side of (4.3). Denote 
this zeroth order approximation by/ (2). Carry out 
the indicated integrations to obtain j x (z) and adjust 
G so that/! (+A)=0. 

2. Repeat, using/i (z) to generate f 2 (z). 'Readjust 
Cso that/2 (+h) = 0. 

In principle this process may be continued indefi- 
nitely with the formula for the iVth approximation 
being 



/i 



-y ff _ 1 («)w(a,r)]«if V 



(4.4) 



subject to /jv(+A)=0, which defines C N , the Nth 
approximation to C. 

Equation (4.4) is formally true for any w(z, f). 
However, it will clearly be to the advantage of the 
investigator to make a choice which results in rapid 
convergence of the iterative process. This matter 
is the raison d'etre of much of the literature of 



linear antenna theory. A choice which leads to 
manageable integrals is 



w H (z,t)=f 



Then 



[a 2 +(z-m m 
C +h d£ 

fe(a)= J-»[0»+(«-f)T*" 



(4.5) 



(4.6) 



The kernel g(z-£) as given by (2.14) and (2.15) 
is difficult to handle. A manageable but crude 
approximation is 



0(z-f)=- 



-s*\'-t\ 



■t\ 



(4.7) 



If g(z-{) is approximated to this order it is ap- 
propriate to approximate w H (z-£) to the same order 
when it is used in the integrals on the right hand 
side of (4.4), but not in the calculation of ^ H {z). 
Then 



In{z)- 



1 



lfe(*) 



{ p <*>-f,- 



(fla-H-fl-/^. 



l*-ri 



=&>*} 



(4.8) 



becomes the fundamental equation of the iterative 
program. Equation (4.8) has been used by Hallen 
in the investigation of linear antenna theory [6]. 

Equation (4.8) can be criticized on two counts. 
First, the approximation involved in replacing R 
by simpy \z-£\ is quite severe. Secondly, w(z,£) 
was chosen for its simplicity rather than according 
to the requirement that the iterative program 
produce good results in low order. King and 
Middleton improved the iterative procedure out- 
lined above in that they used a kernel distance 



fii=[a»+(«-r)»J 



(4.9) 



This is a better approximation to R than is the 
quantity \z-£\. 

The second modification introduced by King and 
Middleton comes from considering the combination 
of integrals on the right hand side of (4.3). Denote 
these by 

<?(s)= P [/(r)y(8-r)-/(*)w(*,r)]#- (4.10) 

J -h 

Q(z) can be rewritten as 

«(*)= P \f(t)-MW(z,!)]g(z-t)d!> (4.11) 

J -h 

A W (z, f) which makes the integrand of (4.11) 
vanish is 



w (z, n=/ (r)//(*). 



(4.12) 
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Of course, one cannot know the desired W(z, f) 
because /(g) is not yet known. One can, however, 
approximate W(z, f) by making use of a fair low 
older approximation to f(z). A suitable choice of 
f(z) in this case is the sinusoidal approximation to 
antenna current. 3 

Instead of W (z, f) King and Middleton introduce 



W (o .x _ sin Ar(A-lfl) 



with the attendant 



Jh p -jhR, 



fli 



(4.13) 



(4.14) 



Using these definitions the interative program is 
based on 



w-i&iw-L 



ikR, -\ 



lU-itt)-fN-i(z)W K (z, w-p-di y. (4.i5) 



Ri 






Equations (4.8) and (4.15) are not quite the forms 
used by Hallen and King-Mid die ton in their compu- 
tational programs. To appreciate the need for some 
improvement consider the boundary condition 
j( + h)=0 applied to (4.4) from which (4.8) and 
(4.15) were developed. The result of applying the 
boundary condition is 



Q=P(h)-f H J N . 



i(f)ff(A-f)rff- 



(4.16) 



(The term involving /y_i (z) drops out because 
the boundary condition is applied at each iteration.) 
Now, P(h) contains C cos kh and if kh=wl2, the 
constant C disappears completely and (4.16) cannot 
be satisfied. A revision of the theory to overcome 
this defect can be made by subtracting (4.16) from 
(4.4) to obtain 

^(^)-^{p(^-^w+£y^-i(f)^(/i-f)^ 

-§ h U N -^)g{z-$)-fN-i{z)w(z, f)]rff|- (4.17) 

The condition f N (+A) = is always automatically 
satisfied by (4.17). Equation (4.16) is forced to hold 
for all N so that C is determined by 

0=P(h)- f fM)<Ah-0<n- (4.18) 

J -h 

A zeroth order approximation which satisfies 
boundary conditions can be obtained directly from 



3 A choice of/(z) to serve in constructing W(z, f) and tf-(z) need not influence the 
choice offo(z) which is used to start the interative process. 



(4.17) by omitting the integrals. Then 
P(z)-P(h) 



/•(*) = 



Kz) 



(4.19) 



Additional notation and definitions may be in- 
vented so that the result of the iterative process can 
be cast into series form. The series can be designed 
to lead off with f (z) as given by (4.19). A suffi- 
cient basis for the remaining part of our discussion 
has been displayed at this point. 

A noteworthy criticism of the use of R x has been 
present by Gans [7]. Gans correctly points out that 



r. 



m) 



-ikRi 



Ri 



■ df=(7eos kz+~ sin k\z\ (4.20) 



is not a true equation because the right hand side 
has discontinuous derivatives in z whereas the 
derivatives of the left hand side are continuous. 

Hallen claims immunity from Gans' criticism on 
the grounds that he uses the (list mice \z— f|. It 
seems to us that this practice raises another diffi- 
culty. A development similar to that used in exam- 
ining the singularity in /(f) shows that g{z — $) is 
logarithmic near z— f=0. The approximate kernel 
\z— f| -1 exp (—jk\z—£\) lias an entirely different kind 
of singularity. It is extremely doubtful if the 1 In I lei i 
theory develops a solution to the original integral 
equation. 

These considerations lead to a definite statement 
that 



f(z)^Limf N (z) 



(4.21) 



where f N (z) is taken from either the Hallen or King 
Middleton form of the theory and/(s) is the correct 
solution of (4.3). Even though one is compelled to 
this conclusion, it is completely irrelevant because 
it is entirely possible for 



/OO^/irOO 



(4.22) 



in low order. If any confusion exists about this 
matter it is because many writers (including some 
authors of senior and graduate level texts) begin 
their discussion directly with eq (4.20). The only 
correct procedure is to formulate a problem which is 
soluble in principle and to introduce judicious 
approximations as needed during the course of 
solution. 

Iterative solutions have the disadvantage t hat they 
become tedious in second and third order even if 
approximations are made. High-order solutions with 
approximate kernels are not even desirable as they 
may have nothing to do with the original problem. 
When all is considered, one needs to know over what 
range of antenna parameters such solutions can be 
used with confidence. Obviously, the kernels used 
in Hallen and King-Middleton theory approximate 
g(z—£) over a range comparable to antenna length 
only if h/a is large. The point at which these theories 
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break down is somewhat arbitrary since it must 
depend upon an arbitrarily selected amount of tol- 
erable error. 

5. Fourier Series Solutions of Storm and 
Zuhrt 

One disadvantage of the iterative solutions dis- 
cussed in section 1 is that they become extremely 
tedious beyond second or third order. Storm at- 
tempted to invent a theory which could be extended 
to higher order [8]. Unfortunately, he introduced 
the kernel distance R 1 =[a 2 +(z — f) 2 ] 1/2 and studied 



/: 



my 



-jkRi 



-df=(7cos kz-\-~ sin k\z\ 
hi & 



(5.1) 



which has, in fact, no solution at all. "Solutions" 
to (5.1) are physically meaningful only if they apply 
to sufficiently thin antennas and if they are re- 
stricted to low order. In addition, Storm's theory 
contains fundamental errors which invalidate his 
solution no matter what kernel distance is used. 
Storm expands the unknown function in the form 

/(f)=J?sin&(/Hrl)+ &n cos (2»+l)xf/2fc. (5.2) 

n = 

Since /(f) is an even function the expansion also 
represents /(f) in the range— A < f <0. 

Following Storm, we seek to determine the co- 
efficients B and F n of the above expansion. The 
result of substituting (5.2) into (5.1) is 

M (z)B+ XI S n {z)F n =C cos fe+i sin k\z\ } (5.3) 



Table 1 



n= 



where 



sinKA-lrD^gptt, (5.4) 

S n (z)= P 

J -h 



(2n+l)7rf<r M i , e E . 

cos- — kr 1 - —Fr-d£. (5.5) 



Bi 



The integrals in (5.4) and (5.5) are somewhat 
difficult to evaluate unless approximations are made. 
Storm replaces the kernel with \z — f | _1 exp (—jk\z — 
f|) outside the range z — 5a<f<2 + 5a. Inside this 
range he replaces the kernel by [a 2 +(z — f) 2 ]" 1/2 . 
These ranges are ambiguous if z is within 5a of the 
ends of the antenna since f must also be restricted 
to the range —h<£<h. Presumably we are not to 
consider values of z too close to the ends of the 
antenna in what follows. Once the approximations 
are made, the integrations required in (5.4) and (5.5) 
can be performed in terms of elementary functions. 
We shall omit the details and return to considera- 
tion of (5.3). 

Storm explicitly satisfies (5.3) at JV+2 points and 
obtains iV+2 equations in N+2 unknowns which 



kh 


Storm: 5 point 
calculation 


King-Middleton 


tj-,2 

7T 


Z* = 81.54-;44. 5 
Z* = 1162 -j 1354 


Zjcjf=81.5+j43.4 
Z K m= 1000 -j 1350 



are B, C, and N of the F n . He has performed cal- 
culations with N=0, 1, 2, and 3 for both full and 
half wavelength antennas with h/a =904. The 
agreement with King-Middleton iterative theory is 
remarkable. A comparison between the impedances 
from the latter theory with those from Storm's 
five-point calculation is made in table 1. King- 
Middleton data used in this comparison are second 
order except for the resistance in the half wavelength 
case which is third order. 

In spite of the success of Storm's calculation, his 
theory breaks down in higher order or for smaller 
h/a ratios. The difficulties are made more evident 
if we specify the unknown function by the expansion 

/(M& cos &±P«+%F. cos-^ 1 ^ 



2h 



2h 



(5.6) 



where the A n are the first N terms of an expansion 
of sin k(h— |f|). As such the A n are known ex- 
plicitly. It is clear that (5.6) approaches (5.2) as 
N becomes large. If (5.6) is used for /(f), (5.4) must 
be changed to 

(6.7) 

The definition of S n (z) does not change. The new 
equation to be satisfied is 

M'o(z)B+"i£ S n (z)F n =C cos kz+hmk\z\. (5.8) 



Set aside, for the moment, the question of how C 
is to be determined and consider the solution in terms 
of C. A set of linear simultaneous equations for the 
unknowns B, F , F ly . . . F N is obtained by satis- 
fying (5.8) at (N+l) values of z. The matrix of 
coefficients is 



M' ( Zl ) 

M' (z 2 ) 



So(Zi) 

So(2 2 ) 



Si(*k) 



S N (z 2 ) 



M' (z N ) S (z N ) S^Sn) 



S N (z N ) 
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By comparing (5.5) and (5.7), one sees that the 
M'o(z m ) are linear combinations of the S n (z m ). The 
determinant of the above matrix is identically zero 
for any N. If the left hand column is replaced by 
M (z m ) so that all of B sin k(h—\$\) is used, the de- 
terminant of the matrix must approach zero as N 
becomes large. No high order solution can be ob- 
tained unless one of the unknowns is assigned an 
arbitrary value. A likely candidate for the assign- 
ment is B which can be selected from the elementary 
induced emf theory according to B=(jZ l47r)I m . 

Evidently, Storm did not become aware of these 
difficulties because he included C among the un- 
knowns. If this is done and the matrix is augmented 
by appropriate bordering elements, a system of 
equations with an unique solution for every N results. 
However, the value of C which is obtained may not 
be the correct one. Storm's expansion for f(z) is 
identically zero at the boundaries. If the coefficients, 
F n , of thai expansion are calculated in terms of a 
spurious C it may happen that the expansion (for 
large N) does not approach zero at the boundaries 
even though it is identically zero at z=±h. The 
only way out of this dilemma is to expand the current 
in functions which are not identically zero at z=±.h 
and then use the boundary condition /( ± A) = to 
determine C. 

An actual calculation to support these 4 criticisms is 
worthwhile. In performing the calculation we used 
the theory of the next section which is equivalent to 
a corrected version of Storm's theory even though 
several of the technical details are quite different. 
Storm was followed to the extent that /(f) was ex- 
panded in the set {cos (2w+l)7rf/2A}, and C was 
treated as an independent unknown. A structure 
with h/a=60 was selected for study. The results 
are shown in figure 2 which displays the real and 
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Figure 2. Impedance versus order of solution (Storm's pro- 
cedure) . 

imaginary parts of antenna impedance as a function 
of the order of solution. Corresponding values of C 
are shown in figure 3. The result of mistreating the 
parameter C is easily observed from these curves. 
We have conjectured that the amplitude of oscilla- 
tion in the results may decrease with increasing h/a. 
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Figure 3. The parameter C versus order of solution (Storm's 

procedure) . 

If this is so it explains why Storm obtained good 
results with a defective theory. 

Zuhrt considered the problem of developing f(z) 
in a Fourier series from a somewhat different point 
of view [9]. Instead of considering a single isolated 
antenna, he formulated a boundary value problem 
for an infinite collinear array of such antennas 
spaced along the 2-axis with centers at z= ± nd, where 
d^>h. Each unit is center-driven by a potential 
V n =(—l) n V impressed across a gap of finite width. 
Ultimately he allows d to approach infinity and the 
gap width to become small. In this limit only the 
center unit remains and his theory represents the 
simple tubular model. This approach is unneces- 
sarily intricate. Zuhrt's final equation can be de- 
rived directly from the Hallen integral equation. 
Equation (3.3) may be written 



»<« 



■""-if. 



J„ (/3a) H? (pa) cos a(z 



■fida. 

(5.9) 



The term cos a(z— f) in the integrand can be ex- 
panded and the term sin a z sin af omitted since 
/(f) is even for a center-fed antenna. Therefore, we 
can consider 



/(f) K(z,£)d?=C cos kz+- sin k\z\, 



(5.10) 



where 



/» 00 

K{s,i) = —j J J (Pa)Hi?> (Pa) cos az cos afdf. 

(5.11) 

Applying the operator L z =b 2 /dz 2 -{-k 2 to both sides 
of (5.10), one obtains 

P f{S)L z K{zJ)dS=kh(z). (5.12) 

J -k 
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Now expand /(f) as 



/(fHsV.cosfinf, 

w=0 



where 



(2n+l)x 



2& 



=H n . 



(5.13) 



(5.14) 



Substituting the assumed current expansion into 
(5.12), one obtains 



n=N C rh 
n=0 V^J-7 



cos H n ?L z K(z, fid? 



\ F «= 



k8(z). (5.15) 



Now multiply both sides of (5.15) by cos H p z, and 
integrate on z over the range —h<z<h to obtain 



•g- ri r» 



cos H p zdz 



j h ^osH n {L 2 K^,0d{\F n = l. (5.16) 

Equation (5.16) generates an infinite set of linear 
simultaneous equations as the indexing parameter, 
#, is allowed to range from zero to infinity. The 
coefficients are given by 

Z pn =-r cos H p zdzl cos H n $L z K(z,$)d$. (5.17) 

Further reduction is accomplished by substituting 
(5.11) for K(z, f) and carrying out the indicated 
operations. When this is done, 

z Pn = -±j(-iy +n H p H n 

J o mHP(ftaW-c?) cos 2 ak f . 



/; 



This is the same as Zuhrt 's formula except for 
trivial differences in notation. 

Zuhrt obtains an Nth order theory by truncating 
the infinite scheme at Nth order. The integral 
which defines the matrix elements is difficult to 
evaluate. Zuhrt resolves this difficulty by resorting 
to graphical integration. Each coefficient has a real 
and an imaginary part, so that an iVth order theory 
requires 2N 2 graphical integrations, a formidable 
amount of labor even for small N. 

6. Further Development of Fourier Series 
Solutions 

The kernel of HalleVs integral equation, 

-jkR 



so?- r)= 



where 



2xJ_, 



B 



d4» 



(6.1) 



i2 = [4a 2 sin 2 0/2+(2-f) a p, (6.2) 

represents the radiation at any value of z on a 



cylinder of radius a from a ring source of radius 
a located at f. As such it is meaningful in the 
entire domain — oo< z<C + oo . However, the integral 
equation is valid only on the range — A<£<+A and 
operations on the source coordinate, f , are restricted 
to the same range. Hence, a special expansion of 
g(z — £) for — h<z^<-\-h is desirable. A frontal 
attack on the problem of obtaining such an expansion 
has been made by Bohn in an investigation of a 
theoretical model which is quite different from the 
one used here [10]. The point of interest at present 
is his method of handling g{z — $) which is the kernel 
of an integral equation occurring in his theory. 
In our notation Bonn's expansion is 



ff(s-f) = Z! y £,G nm e"***i h - m *W 



(6.3) 



for — A<2,f<+A, where the coefficients of the 
Fourier series are given by 



G nm =^ 2 j\J h h g^-t)e- j ^ z/h - m ^ /h) dzd^ 
Substituting (5.9) for g(z— £), (6.4) becomes 

J + 00 
- 00 

Jo(t3a)Hk 2) (Pa) sin (mr—ah) sin (rmr—ah)da. 
(nir—ah)(mir—ah) 



(6.4) 



(6.5) 



This expression exhibits the difficulties involved 
in a direct attack. The reader will appreciate that 
the integration is not trivial. Bohn evaluates (6.5) 
approximately by means of ingenious distortions of 
the contour of integration. Details will be omitted 
here, it being sufficient for our purpose to note that 
no investigator to date has been able to obtain the 
exact coefficients of a double Fourier series repre- 
sentation of g(z — f). 

Fortunately, it is possible to avoid the double 
Fourier series representation entirely. It is only 
necessary to make use of the fact that z and f enter 
only as the square of their difference. Thus, g(z — £) 
is not a general function of (z, f). To proceed, let 



{=(«-«■). 



(6.6) 



If z and f are separately in the range ~h<z^<h, 
then £ is in the range —2h<£<2h. Since g{£) 
contains only £ 2 , a cosine series in the range to 
2h will suffice. Thus, we seek an expansion of the 
form 

T\ m= o o 

fftt ) =^r+ 2 D m cos mrtfth, for < £< 2h. (6.7) 

^ ra = l 

The fundamental formula for any coefficient is 

i r 2h 
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Direct integration of (6.8) is difficult. An indirect 
method can be constructed by writing 



0te)=0i(€)+02($), 



(6.9) 



where gi(£)=g(£) in the range 0<£<2A, and is zero 
outside this range, g 2 {£) is equal to zero in the range 
0<£<2A, and is identical to g(J£) in the range 
2h<£<oo. Transposition of (6.9) gives 

fttt)=fifft)-to(€). (6.10) 

The symbolic cosine transform of (6.10) is 

6 l (a)=6(a)-6 2 (a). (6.11) 

The transform of g x (£) and its inverse are given by 



and 






g x (£) cos ajrff, 



(a) COS afyla. 



(6.12) 
(6.13) 



Comparison of (6.12) and (6.8) shows that the 
coefficients of the cosine series expansion for #i(£) 
are simply proportional to sample values of its 
cosine transform, 6i(a). Thus, 

6(a) can be immediately identified from (5.9) as 

G(a)= -]^lJo(Pa)B?> (/3a). (6.15) 

Incidentally, an asymptotic expansion of G(a) shows 
that, for small £, g(£)~ln £. Since the singularity 
in #(£) is no worse than logarithmic, the integrals 
defining the D m exist. 

Determination of G 2 (a) is tedious but not difficult. 
The definition of 6 2 (a) is 



G 2 (*) = 



VtL 



g (£) cos a##, 



(6.16) 



which can be expanded in a highly convergent 
infinite series. If, for convenience, we set 



and 

(6.1) becomes 



2/ 2 =4a 2 sin 2 0/2 
u 2 =yW, 






(6.17) 
((5.18) 

(6.19) 



When appropriate expansions of the integrand are 
made and the ^-integration is performed, there 



results 



(l(?>y 2 )=~[l- jj y-j-. • •} (6.20) 



The above expression is to be used when £>2h. 
The parameter k is of the same order as 1/h. There- 
fore the terms retained in the square bracket are of 
order (a/A) 2 . The first omitted term is of order 
(a/hy. We shall now omit all of (6.20) except the 
first term. It is not difficult to restore the small 
correction terms later. 
With this omission, (6.16) becomes 



fl'.O 



wu; 



e '**cos a£ 



it 



(6.21) 



Sine and cosine integral functions are defined by 



Ci(x)= 



(6.22) 



COS ./' 



X 



dx. 



With the above, a few trigonometric identities, and 
a few elementary changes of variable, the integration 
of (6.21) follows almost by definition, it being only 
necessary to exercise a little care depending on 
whether «<& or a^>Jc. The case of a=k will require 
special attention. We have 



G 2 (a<k)c 



^lm(k+) 



cQh]-yCi[(k-a)2h] 



+JT-jSi[(k+a)2h\-jSi[(k-a)2h]}- 

(6.23) 

G 2 ( a >k)^-^{Ci[(a+k)2h} + Ci[(a-k)2h} 

-jSi[(a+k)2h]+jSi[(a-k)2h]}- (6.24) 

It is convenient to rewrite 6(a) from (6.15) separ- 
ately for the cases a<Ck and a^>k. Let \P\ = b. The 
phase of jS has been defined so that fi—b if a<Ck and 
/3=-j6 if a>k. Then 

G(a<k)=-j -yjl J (ba)H^(ba), (6.25) 



and 



6(a>k) 



=y - Jo(&a 



)^o(^)- 



(6.26) 



Now ff 2 (a<&), 6 2 (a>k), G(a<k), G(a>k) are 
all singular at a=k. In the first two functions the 
singularity comes from the cosine integral function, 
in the last two Ho m (ba) and K (ba) become singular. 
Since these functions are to be sampled at the points 
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mir/2h, the singularities apparently give trouble if 
A=mX/4. Actually we are concerned only with the 
difference G(a) — G 2 (a) and this turns out to be finite 
and independent of whether a approaches k from 
above or below. The special formula required for 
the a=k case is found from combining (6.23), (6.25) 
and making use of small argument formulas for the 
various functions involved. The latter are tabulated 
for the reader's convenience. 
For small x, 



J (x)~l, Si(x)~0, 

J w yx 



Ci (z)~ln y x. (6.27) 



The logarithmic singularities subtract off in G(a) — 
G 2 (a) and one obtains 



Oi(a=k)c 



Mm 

J2t\ 



4h 



yka 2 



-Ci (4kh)-jSi(4kh) 



y 



(6.28) 



The same result can be obtained by using (6.24), 
(6.26), and the small argument formulas for Io(ba) 
and K Q (ba). 

The degree of approximation in the above formulas 
may be improved by calculating the cosine trans- 
forms of the correction terms in (6.20). The next 
term to be included is 



T 1 (a) = -jka< 



/2 f °° e^L 

YttL ? 



cos ajrff, (6.29) 



which reduces to trigonometric functions, sine in- 
tegral functions and cosine integral functions. All 
higher order correction terms may be similarly 
treated. 

Thus, an expansion of g(z—£) in the form (6.7) 
can be achieved and the coefficients can be calculated 
with any desired degree of accuracy. It will be con- 
venient in what follows to re-define D so that the 
leading term can be included under the summation 
sign. If this is done 



g(z—£) = X) D m [cos mTz/2h cos imrffih 

m=0 

+sin 7mrz/2h sin mw£/2h]. (6.30) 

The sine terms are not needed in the treatment of a 
center-fed antenna. We are then led to consider 



/: 



f(t)K(s,t) dt=C cob kz+hink\z\, (6.31) 



where 



K(z,£) = S D m cos m w z/2h cos rmr f/2A, (6.32) 

w = 

subject to/(±A) = 0. 



We now expand /(f) as 



/(f) = S F n cos riTt/h 



n=0 



(6.33) 



and substitute into (6.31). After the f integration, 



n= oo m= c o 

2—1 2—1 " n^rn) 'nn 
w=0 m = 



cos nnrz/2h=C cos kz-{-= sin k\z\, (6.34) 

7nm=2 cos n 7T £jh cos m t f /2 h d£. (6.35) 

An infinite set of linear simultaneous equations is 
obtained by multiplying both sides of (6.34) by 
{cos (2p + l)Tz/2h} and integrating on z from —h 
to -\-h. The result is 



71 = co W=c o 

S S FnP„ynnfi,m=Cr p + V p , (6.36) 

7i = 771 = 



where 

ft 



'/. 



2 cos [(2^+1)tt0/2^] cos [mirz/2h\dz, 



and 



r p =2 cos kz cos [(2p+l)7r2/2A]<fe, 

^=2 sin kz cos [(2p+l)7r2/2A,]d2. 
Equation (6.36) can be written 



where 



2—1 *■ pn" n ^'p i Vpy 

71 = 



*■ pn 2—1 ^iriinmPp 



(6.37) 

(6.38) 

(6.39) 

(6.40) 
(6.41) 



If (6.41) were actually an infinite sum, many terms 
would be required to satisfactorily approximate each 
T pn and the theory would be laborious except in low 
order. However, (6.41) contains only two nonzero 
terms! To appreciate this fact consider the set of 
functions {cos(mirz/2h)} which appear in the expan- 
sion of K(z, f). These functions are complete on 
0< z<2h. They appeared in the theory because we 
expanded a function of (z— f). Clearly since m is 
either even or odd, this basic set of functions can be 
divided into two subsets { cos rrnrz/h} and {cos (2m + 1) 
Tz/2h}. Both of the subsets are complete and 
orthogonal on the range where they are actually 
used. In a sense we can refer to (6.32) as an over- 
complete expansion of K(z,£). Therefore, the sum- 
mation in (6.41) can be broken into two parts, one 
for m even, the other for m odd, and reduced to 

^pn = h[Pp,2?iD 2 n + yn,2p+lD 2 p+l]' (6.42) 
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The set of equations generated by (6.40) as p 
ranges from zero to infinity can be cast in matrix 
form as 



TF=Cr+v, 



(6.43) 



where T is a matrix of the T pn ; F, r, and v are column 
vectors. 

The method of obtaining simultaneous equations 
from (6.34) used here is formally equivalent to any 
other method which might be used. Consider mul- 
tiplying (6.34) by X v (z) and integrating on z from 

V. ~h to h, where X v (z) is arbitrary. Now X p (z) can 
be expanded in the set {cos {2p-\-\)Tzl2h} so that 
the set of equations obtained by using X v (z) is a 

[ linear combination of the equations represented by 
(6.43). An iVth order colocation scheme is equiva- 
lent to choosing X v (z) from a set of A 7 " delta functions, 
8(z-z p ) withp=l,2,3, . . . . , N. 

The formal solution of (6.43) in terms of F _1 is 
simply 



F=CT- 1 r+T- 1 v. 



(6.44) 



We have not been able to discover a general form for 
F _1 . Consequently, it has been necessary to trun- 
cate the system of equations to finite order and invert 
finite matrices using digital computer methods. 4 
This procedure is a cause of some concern in that 
there are apparently no mathematical theorems 
which justify the assumption that a sequence of 
such inverses converges to the inverse of the infinite 
matrix. However, there is evidence that the pro- 
cedure being used here produces the correct solution 
of Hallen's integral equation. First of all, we have 
examined sequences of finite inverses up to 25th order 
for h/a=60 and kh = T/2. These are well-behaved 
and stable. For a rough definition of stability, we 
shall say that a stable iVth order inverse has been 
found if the elements of an Nih order matrix formed 
by truncating the inverse matrix of an (N-\-M)th 
order solution do not change appreciably as M is 
increased. M will be referred to as the stability 
margin. For the problem at hand our work indi- 
cates that M^3. 

The results obtained for finite antennas are 
thoroughly reasonable when compared to results 
obtained by other methods. From a pragmatic 
point of view there seems to be sufficient evidence 
that the matrix inversion procedure does indeed 
produce a finite number of terms of the correct 
solution to the original integral equation. 

If the elements of the inverse of an JVth order 
truncation of T are designated as H% p the numerical 
solution is 



^(0)=zf#£>(CV p +*g. (6.45) 



The Nth approximation to f(z) given by the 
computational program is 



M*)=T,Fu(C) cos (mcz/h). 

71 = 



(6.46) 



The leading terms of (6.46) are not a good approxi- 
mation to antenna current. Consequently, this 
representation is slowly convergent. To improve 
convergence a good low order approximation io f(z) 
was chosen, and expanded in a series 



x ( z ) ==: ^2 X n cos (nwz/h). 

71=0 



(6.47) 



4 Numerical inversion of the required matrices was accomplished on the IBM 
650, using a library program furnished by IBM Corporation. 



Candidates for the role of x(z) are the classical 
sinusoidal distribution, the zeroth order approxi- 
mation from an iterative solution, the King-Middle- 
ton modified zeroth order approximation, or a new 
low order approximation by K. W. P. King referred 
to by its author as the "quasi-zeroth order approxi- 
mation" [11]. We have used the latter. Once the 
choice has been made and X n have been calculated, 
(6.46) is modified to read 

/»(*)=*(*)+ E [F n (C)-X n ] cos (nirz/h). 

71 = 

(6.48) 

When the boundary condition / (h) = is imposed, 
one obtains an iVth order approximation to (J from 

n=N 

H[F n (C)-X n ](-iy=0. (6.49) 

n = 



By letting N=0, 1, 2, 
generate a sequence. 



., N, (6.49) can be used to 



Go, Gi, 62, . . ., Cn 



(6.50) 



Now an JVth order solution should be regarded as 
the first AT terms of a solution of infinite order, only 
the latter solution involves the Fourier series for 
f(z). The Fourier series is unique and logically it 
must be in terms of C m . The most direct method 
of determining C m is to plot the above sequence 
versus 1/n and extrapolate. Unfortunately, the 
sequence of C values obtained from the boundary 
condition oscillates and the extrapolation is subject 
to large error. This difficulty can be overcome by 
applying a Ceasaro transformation to the sequence 
of C values to form a new sequence which converges 
to the same limit. The transformed sequence plots 
a smooth curve against 1/n, the extrapolated limit 
of which is taken to be C m . A graph of a typical 
treatment of Cis shown in figure 4. 

Our final expression iov j(z) is 

U(z)=x(?)+S [F n (CJ-X n ] cos (mcz/h)- 

n=0 

(6.51) 



581 



.250 



.200 -"~~~"— -- ^__ 




Determination of parameter C. 



This expression does not quite satisfy the boundary 
condition /(+/&) =0. However, except for trunca- 
tion errors which affect only the last two or three of 
the F n appreciably, it represents an estimate of the 
first N terms of a solution of infinite order. Readers 
who may prefer an expression which satisfies the 
boundary condition exactly in JVth order should be 
reminded once again that a finite number of terms 
of a Fourier series provides a least squares best fit 
over the entire range of the function being repre- 
sented. This type of fit is to be preferred over one 
which is identically equal to the function at one 
point. 

The admittance of the antenna is now simply 



n = N 






-*«]}- 



(6.52) 



The function /(f) is proportional to the sum of the 
currents on both the inner and outer surfaces of the 
tubular conductor. It was pointed out earlier that 
only the current on the outer surface and not the 
tube current is to be associated with the experimen- 
tal admittance of the antenna. 

The generic expression for the vector potential at 
a field point (p, <£, z) when p =a is 



where g(p, z; a, f) is given by either 



with 



or 



1 C* 

g(p ) z;a,£)=— J 



-jkR 



R 



d<t> 



R=[ P 2 +a 2 -2pa cos 0+O-f) 2 ] 



1/2 



(6.53) 

(6.54) 
(6.55) 



f* 00 

g( P ,z;a,£) = -j J J o (0a)H^(fi P ) cos a(z-?)da- 

(6.56) 

As a consequence of the ^-symmetry, only the 
^-component of H is different from zero, and is 
given by 



Ht= — dA 2 /bp- 



(6.57) 



The surface current density on the outer surface of 
the tube is obtained from the boundary condition 
nXH=K, and the current distribution along the k 
outer surface is 



where 



P(z)~j h mg\ P ,z-M)d^ (6.58) , 



Sr°( P ,2;a,f)= -2,ra[|^/(p,s;a,r)] ^ • (6.59) , 

If g{p,z;a,£) is expanded in the Fourier series 

fif°(p, 2 ;a,f)= S^S. cos m * { *-S\ 0< \z- f | <A, 

(6.60) < 

the coefficients D^ may be determined by the same 
procedure used in expanding the kernel of Hall en's 
equation. Omitting details, the D^ are given by 
sample values of 



G° l (a)=-jpaJ (pa)H?\Pa) 



f_5_ f " p e~ jkR 



)} 



+■ -i ir I ^~n~ cos a(z—{)d<j>d(z—t) f 

(6.61) 
where R is given by (6.55). 

This leads to the following expression for Fourier 
coefficients of the exterior current: 






iV x{X) cos (mrffthW 

+ 12 [F n -X n ] ( H cos (nwt/h) cos (mrt/2h)d{ V- 

71 = J -h J 

(6.62) 

7. Results 

We have applied the theory to half and full 
wavelength antennas with h/a =60 and 5007r, 
respectively. In each case the calculations were 
performed at 25th order. Correction terms from 
(6.30) were included in the calculation so that each 
matrix element is accurate to one part in 10 5 . 
Graphs of total current and plots of total admit- 
tance versus order are shown (figs. 5 to 12). In 
the latter graphs, an nth order admittance is ob- 
tained by using the first n terms of a 25th order 
solution. 

Comparison of these results with those of the 
King-Middleton theory in terms of impedance are 
shown in table 2. 

King-Middleton impedances in table 2 were ob- 
tained by graphical interpolation of tables given in 
reference [2]. Kesistances of the half wavelength 
structures are from third order solutions. Other 
quantities from King-Middleton theory are second 
order. Our own Z T for the full wavelength antenna 
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with h/a=$0 is obtained by extrapolating the current 
graph to 2—0. The current distribution on a full 
wavelength structure does not have zero slope at the 
origin as do the cosine terms used to describe it. 
Consequently, many cosine terms are required for 
high accuracy, especially in the treatment of a thick 
structure. The inner current correction is entirely 
negligible for structures with h/a=500 tt; it is still 
small for h/a=Q0. 
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Figure 5. Current distribution: kh=7r/2, /i/a = 60. 
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Figure 6. Admittance versus order of solution: kh = ?r/2, 
h/a = 60. 
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Figure 7. Current distribution: kh = 7r, h/a = 60. 
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Figure 8. Admittance versus order of solution: kh = 7r, h/a 

= 60. 
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Figure 9. Current distribution: kh = 7r/2, h/a=500 jr. 
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Figure 10. Admittance versus order of solution: kh = 7r/2' 
h/a = 500 7T. 
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Figure 11. Current distribution: kh = 7r, h/a = 500 w. 
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Figure 12. Admittance versus order of solution: kh=7r, 
h/a = 500 7T. 

It can be seen from table 2 that the amount of 
disagreement between this theory and the King- 
Middleton theory is only about two percent. One 
may expect disagreement of this order or less over the 
entire range of h/a^>60. 







Table 


2 




kh 


hja 


Zt 


z° 


Zkm 


7I-/2 

7T 
7T/2 


60 
60 
500 w 

500 7T 


91.4+j 38.6 
205 -j 382 

79. 7-j 42. 9 
1646 -j 1768 


92.5+;' 40.6 
205 -j 380 


91. 4+j 41. 5 
206 -j 380 
SO. 3-}' 43. 4 






1625 -;1744 







It can be seen from the graphs of admittance versus 
order that, except for thick full wavelength struc- 
tures, 25 terms are excessive for A/a>60. Thus, a 
great deal of margin for the study of thicker antennas 
is inherent in 25th order solutions. Studies of full 
wavelength antennas with hja < 60 will require more 
than 25 terms or further modification of some of the 
technical details of the theory. 
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